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This  paper  discusses  the  numerical  simulations  of  selected  high-speed  aeropropulsive 
applications  including  axisymmetric  supersonic  base  flowfields,  with  and  without  jets,  and 
three-dimensional  supersonic  hot  jets.  Comparisons  are  made  using  a  compressibility 
corrected  k-s  turbulence  model  and  an  Explicit  Algebraic  Stress  Model  (EASM).  The 
present  EASM,  referred  to  as  EASM/J,  has  been  modified  to  correctly  predict  the  mixing  of 
basic  jet  and  free  shear  flows.  A  structured,  parallel,  finite-volume,  Navier  Stokes  code, 
CRAFT  CFD  is  used  as  the  flow  solver.  In  the  present  study,  the  EASM/J  analysis  assumes 
a  fixed  (equilibrium)  value  of  production-to-dissipation  ratio.  The  effect  of  the 
compressibility  correction  as  implemented  in  CRAFT  CFD,  in  predicting  the  flow 
reattachment  location  and  base  pressure  distribution  is  investigated.  For  base  region 
studies,  the  k-s  and  EASM/J  models  predict  comparable  results,  but  base  pressure 
distributions  are  closer  to  the  experimental  data  without  using  the  compressibility 
correction.  For  the  jet  with  base  study,  the  predictions  with  and  without  the  correction 
bound  the  data.  For  the  three-dimensional  cases,  there  is  no  substantive  difference  in  mean 
flow  predictions  using  k-c  and  EASM/J. 


Nomenclature 

ui  =  Cartesian  components  of  velocity 

t  =  time 

p  =  density 

juL  =  laminar  viscosity 

T..  =  shear  stress  tensor 

k  =  turbulent  kinetic  energy 
s  =  turbulent  dissipation  rate 
juT  =  turbulent  eddy  viscosity 

P/s  =  ratio  of  production  of  turbulent  kinetic  energy  to  turbulent  dissipation  rate 
by  =  anisotropy  stress  tensor 

S..  =  strain-rate  tensor 

lJ 

W..  =  rotation  rate  tensor 

a.  =  coefficients  calculated  in  EASM/J 

D  =  the  combined  effect  of  turbulent  transport  and  viscous  diffusion 
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I.  Introduction 

This  paper  compares  the  performance  of  an  aero-propulsive  variant  of  the  EASM  turbulence  model 2’3,  to  that  of 
the  k-c  turbulence  model  4  used  in  our  jet  and  aeropropulsive  studies,  for  a  series  of  test  cases  of  increasing 
complexity.  The  k-s  model  uses  the  original  Launder  et  al.  5  coefficients  and  contains  a  compressibility-correction  4 
found  to  adequately  analyze  a  variety  of  basic  high-speed  jets  and  free  shear  flows.  The  EASM  utilizes  this  k-c 
model  for  obtaining  needed  turbulent  kinetic  energy  and  time-scale  information  and  makes  use  of  a  recalibrated 
pressure-strain  correlation  coefficient  to  perform  comparably  to  the  k-c  model  for  basic  jet  and  free  shear  flow 
mixing  2,3 . 

EASM  models  have  several  potential  advantages  over  k-c  models  that  include  realizability  (in  a  more 
fundamental  framework  than  achieved  via  just  adding  realizability  modifications  to  k-c)  and  a  better  prediction  of 
turbulent  stresses,  in  particular  normal  stress  anisotropy.  As  we  move  from  simple  jets  and  free  shear  flows,  in  this 
paper  we  examine  what  improvements  over  k-c  model  are  afforded  by  EASM  for  a  variety  of  more  complex 
aeropropulsive  flows.  This  paper  deals  with  issues  of  flow  structure  as  exhibited  by  the  improved  behavior  with 
EASM  for  a  missile  at  angle-of-attack  with  cross-flow  separation  (Figure  1).  Issues  of  anisotropy,  as  related  to  jet 
noise  problems,  are  the  subject  of  a  related  paper  at  this  meeting  6.  We  have  used  the  basic  c-equation  in  all  these 
studies,  realizing  that  extended  variants  will  be  needed  for  complex  flows  with  curvature,  swirl,  etc.,  but  reserving 
the  exploration  of  such  corrections  for  a  future  paper.  Both  the  k-c  and  EASM  models  have  low-Re  and  wall 
function  capabilities  to  treat  near-wall  regions,  and  the  overall  equations  will  be  described  in  the  Section  III. 
Specific  problems  that  have  been  analyzed  include  that  of: 

1 .  Base  region  behind  a  cylinder7; 

2.  Propulsive  jet  into  a  moving  stream  with  a  large  base  8; 

3 .  Lateral  slot  j et  inj  ection  9; 

4.  Propulsive  jet  with  lobes  in  the  exhaust;  and, 

5.  Dual  supersonic  interacting  jets. 

The  authors  and  coworkers  have  been  involved  in  the  extension  and  unification  of  turbulence  models  for 
application  to  a  broad  variety  of  aero-propulsive  flows  ranging  from  laboratory  jets  and  aircraft  exhausts,  to  scramjet 
combustors  and  rocket  plumes  10,1 1,12 .  Thermal/species  turbulent  transport  issues  can  be  quite  significant  in 
aeropropulsive  flows,  and  are  also  the  subject  of  a  related  paper  at  this  meeting13. 


(a).  Comparison  of  stagnation  pressure  contours  at 
selected  stations  along  missile  at  a=14-deg 


(b).  Predicted  pitot  pressure  improvements  at  X/D=l  1.5 
with  EASM1. 


Figure  1. 


II.  Background 

Most  turbulence  models  within  a  RANS  framework  are  linear  eddy  viscosity  models,  which  assume  the 
Boussinesq  hypothesis  to  be  valid.  For  some  of  the  complex  turbulent  flows  of  engineering  importance,  this 
assumption  may  fail.  Full  Reynolds  stress  models  are  more  accurate  for  modeling  such  flows  as  they  do  not  assume 
the  turbulent  shear  stresses  to  be  proportional  to  the  rate  of  mean  strain,  and  thus  can  be  expected  to  provide  better 
predictions  for  flows  with  sudden  changes  in  the  mean  strain  rate,  or  with  effects  such  as  streamline  curvature. 
However,  keeping  in  view  the  computational  effort  needed  to  implement  these  models,  a  more  practical  approach,  is 
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that  of  using  Algebraic  Stress  Models  (ASM).  Unlike  full  Reynolds  Stress  Models,  which  require  solving  the  full 
partial  differential  equations  for  each  of  the  stresses,  Algebraic  Stress  Models  provide  a  nonlinear  relationship 
between  the  Reynolds  stresses  and  the  mean  rate  of  strain.  Several  researchers  developed  these  nonlinear 
relationships  in  the  form  of  a  series  expansion  having  the  Boussinesq  assumption  as  the  leading  term.  Early  studies 
by  Rodi  14  deduced  the  nonlinear  algebraic  equation  for  the  Reynolds  stresses  by  simplifying  the  full  Reynolds  stress 
partial  differential  equation  (PDE).  Algebraic  Stress  Models  are  derived  from  the  full  Reynolds  stress  transport 
equation,  and  hence  are  expected  to  perform  better  in  predicting  the  flow  physics.  Initial  attempts  to  solve  ASM 
equations  implicitly  limited  the  allowable  time  step  size  to  obtain  a  stable  solution,  particularly  for  complex  flows. 
Pope  15  developed  a  procedure  for  obtaining  explicit  solution  to  the  proposed  nonlinear  algebraic  equation  for  two- 
dimensional  mean  flows.  Gatski  and  Speziale  16  came  up  with  a  similar  method  for  three-dimensional  flows.  In  the 
present  EASM  model  however,  the  approach  of  Jongen  and  Gatski  17,18  for  two-dimensional  flows  is  used,  and 
applied  for  complex  aeropropulsive  applications.  A  brief  overview  of  the  turbulence  models  used  for  current 
research  is  provided  in  the  following  two  sections. 


III.  Overview  of  Unified  Variant  of  k-s  Model 

The  governing  equations  are  Favre  averaged  (mass  averaged)  Navier-Stokes  equations.  The  closure  to  the  Favre 
averaged  Navier-Stokes  equations  is  achieved  by  relating  the  Reynolds  stresses  to  available  mean  flow  gradients, 
such  as  in  the  Boussinesq  approximation.  This  requires  specification  of  turbulent  eddy  viscosity,  which  presently  is 
determined  using  the  k-c  equations.  The  transport  equation  for  turbulent  kinetic  energy  is  given  by, 
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and  the  transport  equation  for  dissipation  rate  equation  is  given  by, 
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where  juT=C^p —  ,  and  the  coefficients  used  in  the  model  are  C£  =  1.43  ,  C£  =  1.92  ,  CM  =  0.09 ,  ok  -  1.0 ,  and 

cje  =  1.3.  SSk  is  the  modeled  pressure  dilatation  term,  which  accounts  for  reduced  mixing  in  compressible 

turbulent  flows.  The  present  approach  is  a  combination  of  the  Sarkar  model,  and  Zeman  ’  model.  The 
functioning  of  Sarkar19  model  in  a  compressible  flow  is  through  the  dilatational  and  pressure  fluctuation  terms.  The 
original  Sarkar  model  involves  correction  at  all  levels  of  compressibility.  For  lower  speed  compressible  mixing 
layers,  this  produces  unwanted  mixing  reductions  and  smaller  than  experimentally  measured  growth  rates.  The 
Zeman  correction  introduces  a  lag  wherein  the  turbulence  fluctuation  due  to  compressibility  must  reach  a  certain 
threshold  before  affecting  the  dissipation.  However,  the  original  model  of  Zeman  does  not  account  for  the 
additional  effect  of  pressure  dilatation.  The  current  method  combines  the  pressure  dilatation  of  Sarkar  with  the  lag 
concept  of  Zeman  to  reduce  compressibility  effects  for  low  speed  compressible  flows  4.  For  clarity,  a  short 
description  about  the  working  of  each  model  is  provided  here  before  explaining  the  present  method.  The  original 
Sarkar  correction,  which  accounts  for  dilatational  and  pressure  dilatational  effects,  is  given  as 

u[.  +  pu[.  m  -a2MTPh  -  (cq  -a3)M2T  ps  (3) 

where  cq  =  0.5 ,  a2  =  0.4 ,  cq  =  0.2  .  This  correction  was  observed  to  introduce  excessive  dissipation  and  growth 
rate  reduction  at  lower  Mach  numbers.  To  minimize  this  effect  at  low  levels  of  compressibility,  Zeman  introduced  a 
lag  term,  which  eliminates  the  effect  of  compressibility  correction  term  at  low  Mach  numbers. 

Mt  =  max  (Mt  —  A,  0)  (4) 


3 

American  Institute  of  Aeronautics  and  Astronautics 


AIAA-2005-1101 


where  Z  is  an  appropriately  chosen  lag  value  to  eliminate  the  effect  of  compressibility  for  lower  Mach  number  jets. 
Introduction  of  MT  required  recalibrating  ai  coefficients;  and  the  final  compressibility  correction  (c.c.)  in  the 
present  model  is  given  by, 


SSt  =  -a[M;Pk-a1M1rps 


(5) 


where  cq  =  2.5 ,  a2  =  2.0  ,  A  =  0.2  .  The  compressibility  correction  has  been  validated  for  a  range  of  basic  flows 
at  low  and  high  Mach  numbers  4.  However,  some  deficiencies  associated  with  the  use  of  this  correction  were  found 
for  complex  aerpropulsive  flows,  and  results  obtained  with  and  without  using  this  correction  will  be  discussed  in 
Sec.  V.  The  coefficients  used  in  the  turbulence  model  are  applicable  for  high-Re  flows,  and  modifications  are 
necessary  for  improved  prediction  of  the  flow  in  the  vicinity  of  a  wall.  The  next  subsection  discusses  the  near-wall 
model  used  in  the  present  analysis. 

A.  Near-wall  Model 

Several  criteria  need  to  be  satisfied  when  choosing  a  near  wall  extension  such  that  it  is  compatible  with  the  high- 
Re  k-e  model.  In  CRAFT  CFD®,  a  modified  version  of  So,  Sarkar,  Gerodimos,  and  Zhang  (SSGZ)  model  22  is 
employed.  The  damping  coefficients  used  in  the  original  SSGZ  model  have  been  modified  to  blend  with  the  high- 
Re  k-c  model,  and  the  resulting  changes  to  terms  in  the  turbulent  dissipation  rate  transport  equation,  given  by  Eqn. 
(2),  are  provided  here.  The  first  term  in  the  RHS  of  Eqn.  (2),  representing  production  of  dissipation  is  modified  to 
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The  second  term  in  the  RHS  of  Eqn.  (2),  representing  destruction  of  dissipation  is  modified  to 
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Finally,  the  coefficient  of  turbulent  eddy  viscosity  is  modified  to 
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In  the  above  equations,  Rer  =— —  ,  is  the  turbulent  Reynolds  number,  and  Re  =— — y  ,  is  wall-distance  based 

Pls  y 

Reynolds  number.  The  k-s  turbulence  model  discussed  above  assumes  a  Boussinesq  relation  to  relate  the  Reynolds 
stresses  to  the  mean  rate  strain.  EASM/J,  which  uses  a  nonlinear  relationship  between  the  Reynolds  stresses  and  the 
mean  rate  of  strain,  is  discussed  in  the  next  section. 


IV.  A  Brief  Overview  of  EASM  Model 

Derived  from  the  Reynolds  stress  transport  equation,  the  equation  for  Reynolds  stress  anisotropy  tensor  is  given  by, 


Dt 
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The  coefficients  at  are  directly  related  to  the  pressure-strain  correlation  model  by  al=—  —  C2 

2  v  3 


a2=-^(2-C4)  ,a3=-^(2-C3)  ,a4=gz ,  where  C2=0.36,  C3- 1.25,  C4=0.4,  z  is  the  turbulence  time  scale,  k/  s  . 


Obtaining  an  explicit  expression  for  the  anisotropy  tensor  involves  certain  assumptions  to  LHS  of  Eqn.  (12),  and  it  is 
explained  in  the  next  few  lines.  An  implicit  solution  for  the  above  Reynolds  stress  tensor  anisotropy  equation,  given 
in  Eqn.  (12)  is  obtained  by  assuming  weak  equilibrium  condition  on  the  turbulent  stress  anisotropy  first  proposed  by 
Rodi 14.  This  is  expressed  as, 


Dt 


(13) 


and  a  second  assumption  is  made,  which  considers  anisotropy  of  the  turbulence  transport  and  viscous  diffusion  is 
proportional  to  the  anisotropy  of  Reynolds  stresses, 


(14) 


Substituting  Eq.  (13)  and  Eq.  (14)  in  Eq.  (12),  we  obtain  a  locally  implicit  expression  for  the  anisotropy  tensor. 
Following  the  work  of  Jongen  and  Gatski 17,  an  explicit  representation  of  the  anisotropy  tensor  can  be  written  as 
follows. 


b  =  axS  +  a2  (SW  -WS)  +  a3 


The  term  g  ,  which  is  required  for  calculating  coefficient  a4  is  written  as  follows: 


(15) 


(16) 

Speziale  and  Gatski  16  determined  a  fixed  value  for  this  expression,  using  homogeneous  turbulence 
considerations.  Direct  application  of  this  formulation  can  be  problematic  in  the  freestream,  and  requires 
regularization.  Alternately,  the  variable  production  to  dissipation  ratio  (P/s)  model  in  the  EASM  has  been 
developed  further  such  that  non-linearity  is  accounted  for  in  the  calculation  of  coefficient  g  due  to  its  dependence 

on  (P/s).  The  expression  for  g  is  given  as: 
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Further  details  of  the  EASM  formulation  are  provided  in  Ref.  16. 

V.  Computational  Methodology 

The  governing  Navier-Stokes  equations  are  solved  using  CRAFT  CFD®  flow  solver.  CRAFT  CFD®  is  a  three- 
dimensional,  Navier-Stokes,  finite-volume,  structured,  parallel,  Fortran90  code,  with  several  features  available  to 
handle  complex  flows.  For  time  integration,  implicit,  ADI,  up  to  second-order  accuracy  with  Newton  sub-iterations 
is  available.  Also,  LU  factorization  method  and  fourth-order  Runge  Kutta  explicit  scheme  can  be  used  for 
approximating  the  temporal  derivatives.  For  RANS  applications,  a  second-order  accurate  upwind  scheme  is  used  for 
spatial  discretization.  A  Roe  flux  splitting  difference  scheme  is  used  for  the  treatment  of  inviscid  fluxes.  The  grid 
blanking  feature  in  CRAFT  CFD®  facilitates  the  placement  of  boundaries  away  from  regions  with  large  flow 
gradients.  Also,  a  physical  boundary  in  the  flow  field  need  not  be  a  block  boundary,  using  this  grid  blanking 
methodology.  A  domain  decomposition  procedure  is  employed  for  dividing  the  entire  computational  grid  into 
specified  number  of  processors,  and  the  inter-processor  communication  is  facilitated  through  MPI-based 
parallelization  strategy. 


VI.  High  Speed  Aeropropulsive  Applications 

The  results  obtained  in  studying  the  selected  aeropropulsive  flows  mentioned  in  Section  I  is  discussed  here. 
Also,  more  generalized  flow  physics  will  be  emphasized  using  the  EASM  model,  particularly  for  flows  with 
substantive  streamwise  vorticity  where  the  k-e  model  has  known  deficiencies.  This  section  is  divided  into  five  sub¬ 
sections,  each  discussing  the  flow  analysis  using  k-c  ,  and  EASM/J  models  for  a  specific  application. 

A.  Base  Region  Study 

Flow  over  axisymmetric  step  with  a  freestream  Mach  number  2.5  has  been  studied  in  detail.  In  this  section,  we 
discuss  results  obtained  using  both  k-c  and  EASM/J  models  compared  with  the  experimental  data  of  Herrin  and 
Dutton7.  This  particular  case  was  investigated  by  many  researchers  to  gain  a  better  understanding  of  the  complex 
flow  behavior,  which  involves  flow  separation,  shear  layer  development,  flow  reattachment,  recompression  shock 
etc.  A  schematic  of  these  fluid  dynamic  processes  is  shown  in  Fig.  2(a).  The  diameter  of  the  cylinder  (i.e.  step 
height)  is  63.5  mm,  and  the  reader  is  referred  to  Ref.6  for  the  details  of  experimental  set-up. 

The  approach  boundary  layer  profile  1mm  upstream  of  the  step  is  generated  by  simulating  the  upstream 
supersonic  flow  at  freestream  Mach  number  2.5.  At  the  inflow  boundary,  flow  variables  p,  V,  and  T  are  computed 
from  the  experimental  stagnation  conditions,  P0  =  5 15kPa,  and  T0  =  294K.  A  Parabolized  Navier-Stokes  approach  is 
used  for  performing  the  upstream  boundary  layer  simulation.  For  the  axisymmetric  step  case,  no-slip  conditions  are 
prescribed  on  the  step  surface.  Non-reflecting  boundary  condition  is  specified  on  the  outflow  boundary,  where  the 
flow  variables  are  extrapolated  from  the  interior  by  letting  the  first-order  derivatives  of  all  the  flow  variables  normal 
to  the  boundary  to  vanish.  Figure  2(b)  shows  the  single-block  grid  generated,  and  Fig.  2(c)  shows  the  Mach 
contours  obtained  using  k-c  model,  without  compressibility  correction  (c.c.).  Figure  2(d)  shows  the  comparison  of 
the  ratio  of  mean  axial  velocity  to  freestream  velocity  distribution.  The  effect  of  compressibility  correction, 
discussed  in  Section  III,  on  the  flow  features  has  been  studied.  Higher  base  pressure  is  predicted  using  the 
compressibility  correction  as  shown  in  Fig.  2(e),  and  the  flow  doesn’t  turn  down  as  much.  EASM/J  and  k-c  model 
predict  comparable  results  with  and  without  compressibility  correction.  The  functioning  of  compressibility 
correction  is  being  investigated  further,  with  the  possible  need  to  use  modified  forms  of  the  epsilon  equation,  or 
rapid  strain  modifications  to  the  compressibility  correction  term. 

B.  Propulsive  Jet  with  Base  Region  Study 

Considering  a  more  general  configuration  for  study  of  the  base  region  in  a  propulsive  flow,  a  case  of  Petrie  and 
Walker  has  been  studied  by  Papp  et  al.23.  The  experiments  are  performed  for  a  missile  configuration,  with  a 
cylindrical  afterbody  and  a  propulsive  jet  in  a  Mach  1.4  free  stream  flow.  The  flow  conditions  are  P^,  =  42.316  kPa, 
^  =  259. 6K,  P0J  =  2121.516kPa,  T0J  =  361.48K.  The  flow  involves  important  features  like  rapid  expansion  at  the 
base  comer,  occurrence  of  Mach  disk  downstream  of  the  base  etc.  A  grid  adaptive  study  has  been  performed  as  a 
part  of  the  analysis  to  adequately  resolve  the  Mach  disk  behind  the  base.  Figure  3(a)  shows  the  contours  of 
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temperature.  The  base  pressure  distribution  obtained  with  and  without  the  compressibility  correction,  using  the  k-s 
model  are  shown  in  Fig.  3(b).  Similar  to  the  observations  made  for  the  axisymmetric  step  flow  in  Sec  VI-A,  a 
higher  base  pressure  is  predicted  when  the  compressibility  correction  is  used.  With  c.c.,  k-s  and  EASM/J  models 
results  are  compared,  and  they  match  closely.  Our  base  region  results  show  that  when  a  boundary  layer  goes 
through  a  rapid  comer  expansion,  compressibility  effects  are  over  estimated  using  the  c.c.  model  of  Papp  and  Dash4, 
which  was  calibrated  for  balanced  pressure  round  jets  and  simple  shear  layers.  The  use  of  EASM  produces  results 
comparable  to  those  of  k-s  and  does  not  improve  the  simulation.  However,  for  other  complex  flows  such  as  the 
lateral  jet  discussed  below,  the  compressibility  correction  works  well  and  provides  marked  improvements. 


(a)  Schematic  of  supersonic  flow  over 
axisymmetric  step7 


(b)  Single  block  grid  generated  for  simulating 
flow  over  axisymmetric  step. 


(d)  Comparison  of  axial  velocity  distribution  with 
experimental  data  of  Herrin  and  Dutton 1 . 


(c)  Contours  of  Mach  number,  using  k-e  model 
without  compressibility  correction. 


(e)  Comparison  of  of  base  pressure  distribution  with 
experimental  data  of  Herrin  and  Dutton  7. 


Figure  2.  Supersonic  flow  over  axisymmetric  step. 
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(a)  Contours  of  Temperature. 


(b)  Comparison  of  base  pressure  distribution  with 
experimental  data  of  Petrie  and  Walker8. 

Figure  3.  Propulsive  jet  with  base  region  results23. 


C.  Lateral  Slot  Jet  Study 

A  sonic,  underexpanded,  gaseous  jet  from  a 
transverse  slot  issuing  into  a  supersonic  external 
flow  is  the  object  of  this  study,  with  the  data  of 
Spaid  and  Zukoski 9  used  for  comparisons.  The  slot 
jet  stream  consists  of  nitrogen  gas,  and  the  external 
stream  is  air  at  Mach  number  3.5.  The 
underexpanded  jet  penetrates  into  the  external 
stream,  is  turned  by  the  momentum  of  the  external 
stream,  and  it  obstructs  the  external  flow  producing 
a  shock  that  causes  the  boundary  layer  to  separate. 
The  shock/boundary  layer  interaction  upstream  of 
the  step  plays  significant  role  in  this  problem,  and  it 
is  found  that  the  separation  zone  ahead  of  the  jet  is 
strongly  influenced  by  compressibility  effects. 
Performing  a  preliminary  analysis  using  EASM/J, 
the  surface  pressure  distribution  is  compared  with 
the  results  predicted  using  k-e  model  by  Papp  at  al. 
[4]  in  Figure  4.  For  this  analysis,  compressibility 
correction  was  not  used  with  EASM/J,  and  the  trend 


Figure  4.  Comparison  of  base  pressure  distribution 
with  experimental  data  of  Spaid  and  Zukoski9. 


of  the  pressure  distribution  obtained  with  EASM/J  is  comparable  with  the  surface  pressure  distribution  obtained  with 
k-c  model,  without  the  correction.  Using  the  compressibility  correction,  the  size  of  the  upstream  separated  region  is 
substantially  larger  and  compares  well  with  data.  For  a  similar  problem  with  a  helium  jet,  comparisons  obtained 
with  compressibility  correction  were  not  as  good  13 . 


D.  Propulsive  Jet  with  Lobes  Study 

Earlier  studies  of  passive  noise  reduction  concepts  for  military  aircraft  involved  the  use  of  chevrons,  and  micro¬ 
jets  10.  However,  the  amount  of  reduction  in  noise  achieved  was  not  large  relative  to  the  observed  drop  in  the  engine 
performance.  A  recent  and  novel  concept  of  ‘corrugated  nozzle’  has  been  designed  at  NCPA/U.  MISS,  and  the 
numerical  simulations  for  this  corrugated  nozzle  concept  \  are  discussed  in  this  section.  These  corrugations  are 
lobes  installed  into  the  divergent  portion  of  the  nozzle,  and  extend  from  the  throat  to  the  exit  of  the  nozzle.  The  exit 
views  of  a  baseline  nozzle,  and  a  nozzle  installed  with  6-lobes  are  shown  in  Figure  5.  The  presence  of  lobes  in  the 
nozzle,  promotes  strong  streamwise  vorticity  in  the  developing  jet  region,  and  thereby  enhances  mixing  of  the  jet 
flow  with  the  surrounding  freestream  flow.  Numerical  simulation  of  the  six-lobed  nozzle  has  been  performed  using 
k-c  model  turbulence  model.  For  better  prediction  of  the  turbulent  flow  structure,  and  to  assess  the  ability  of 
EASM/J  in  fulfilling  this  objective,  EASM/J  analysis  was  also  performed. 

Subsonic  flow  has  been  prescribed  at  the  inflow  boundary  specifying  total  pressure,  P0,  and  total  temperature,  T0. 
Non-reflecting  boundary  condition  is  specified  on  the  outflow  boundary,  where  the  flow  variables  extrapolated  from 
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the  interior  by  letting  the  first-order  derivatives  of  all  the  flow  variables  normal  to  the  boundary  to  vanish.  The 
viscous  walls  are  treated  using  compressible  wall  function  approach,  available  in  CRAFT  CFD®.  On  the  farfield 
boundaries,  freestream  boundary  conditions  are  specified,  allowing  for  entrainment  of  the  jet.  Assuming  symmetry, 
the  simulation  has  been  performed  for  a  half-corrugation. 

The  baseline  k-s  model  used  can  adequately  address  the  mean  flow  effects  of  turbulent  structures  but  cannot 
portray  the  Reynolds  normal  stress  anisotropy  observed  in  laboratory  jets.  EASM  has  the  ability  to  better  resolve 
turbulent  flow  structures  in  complex  vortical  flows  compared  to  the  basic  k-e  model.  This  is  an  outcome  of  its 
development  to  satisfy  turbulence  realizability  constraints  and  results  in  local  modifications  to  Cp  based  on  the 
mean  flow  strain  and  vorticity  behavior.  This  is  shown  more  clearly  by  the  turbulent  kinetic  energy  contours  along 
each  of  the  jet  symmetry  planes  in  Figures  6  and  7.  In  the  immediate  vicinity  of  the  nozzle  lip,  velocity  strain  values 
are  significant,  and  the  EASM  is  observed  to  predict  lower  turbulent  levels  in  this  region,  which  is  felt  to  be  more 
realistic.  Since  the  developing  region  is  finite  and  affected  by  upstream  conditions,  such  small  modeling  details  can 
impact  downstream  flow  behavior  until  self-similarity  is  achieved. 


tke 

i  ,  ,  m  ■  m — 

0.0e+000  1.4e+003  2.9e+003  4.3e+003  5.7e+003  7.1e+003  8.6e+003  1.0e+004 


Figure  6.  Turbulent  kinetic  energy  contours  showing  reductions  for  6-lobe  nozzle  with  EASM/J  (bottom  half) 

and  without  EASM/J  (top  half):  (a)  in  plane  of  lobes. 


Figure  7.  Turbulent  kinetic  energy  contours  showing  reductions  for  6-lobe  nozzle  with  EASM/J  (bottom 
half)  and  without  EASM/J  (top  half):(b)  symmetry  plane  between  lobes. 
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E.  Dual-Jet  Interaction  Study 

In  the  present  study,  flow  through  a  twin-jet  nozzle  at  overexpanded  flow  conditions  has  been  simulated  to  fulfill 
two  objectives.  Firstly,  to  study  the  effect  of  plenum  design  on  the  quality  of  flow  at  the  nozzle  exit.  Secondly,  to 
study  the  flow  features  of  jet  alone.  Following  the  first  phase,  we  observed  that  the  flow  is  quite  non-uniform  as  it 
reaches  the  nozzle  exit,  and  the  plenum  design  was  modified  to  provide  a  much  more  uniform  flow.  Only  the  details 
of  second  phase  of  the  simulation  for  nearly  uniform  flow  nozzle  exit  plane  conditions  are  discussed  in  the  paper. 

Subsonic  flow  has  been  prescribed  at  the  inflow  boundary  specifying  total  pressure,  P0,  which  corresponds  to 
NPR  (Nozzle  Pressure  Ratio)  =  3.94;  and  total  temperature,  T0  =  1022K.  Non-reflecting  boundary  condition  is 
specified  on  the  outflow  boundary,  and  viscous  walls  are  treated  using  compressible  wall  function  approach.  On  the 
farfield  boundaries,  freestream  boundary  conditions  are  specified,  allowing  for  entrainment  of  the  jet. 

The  key  features  of  the  supersonic  jet  issuing  from  the  nozzle  into  still  air  are  shown  in  Mach  number  contours 
of  Figure  8.  The  jet  is  overexpanded,  and  a  Mach  disk  occurs  at  approximately  one  diameter  downstream,  across 
which  the  supersonic  flow  became  subsonic.  The  shear  layers  emanating  from  the  lip  of  the  nozzle  merge  at 
approximately  20  diameters  downstream,  which  forms  the  end  of  potential  core.  This  region  is  marked  by  a  peak  in 
the  turbulent  kinetic  energy  levels,  as  observed  in  the  centerline  variation  of  turbulent  kinetic  energy  in  Figure  9.  To 
get  a  better  understanding  of  the  turbulence  structure,  the  Reynolds  normal  stresses  are  plotted  in  Figures  10-12. 
Major  part  of  the  energy  exchange  from  the  mean  flow  to  the  turbulence  is  observed  to  occur  through  the  streamwise 
Reynolds  normal  stress.  Unlike  the  isotropic  Reynolds  normal  stresses  predicted  using  k-s  model,  there  are  large 
differences  in  the  values  of  the  normal  stresses  predicted  with  EASM/J.  This  is  a  characteristic  feature  of  EASM/J 
as  it  accounts  for  non-isotropy  in  the  evolution  of  stresses  via  non-linear  polynomial  tensor  representation  of 
Reynolds  stresses.  The  EASM/J  calculations  are  felt  to  be  more  realistic  and  results  will  be  validated  by  upcoming 
experiments  at  FSU  by  Krothapalli  and  coworkers. 


Figure  8.  Mach  number  contours,  showing  the  flow  structure  predicted  using  k-s  model. 
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Figure  9.  Turbulent  kinetic  energy  contours  obtained  using 
(a)  k-s  model  (b)  EASM/J  model. 


(b)  EASM/J 


Figure  10.  Streamwise  Reynolds  normal  stress,  u'u'  contours  obtained  using 
(a)  k-s  model  (b)  EASM/J  model. 
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(b)  EASM/J 


Figure  11.  Transverse  Reynolds  normal  stress  contours,  vV  obtained  using 
(a)  k-s  model  (b)  EASM/J  model. 


(b)  EASM/J 


Figure  12.  Spanwise  Reynolds  normal  stress  contours,  w'w'  obtained  using 
(a)  k-s  model  (b)  EASM/J  model. 
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VII.  Conclusions 

The  application  of  an  explicit  algebraic  stress  model  and  its  usage  with  a  turbulence  compressibility  correction 
for  simulating  complex  aeropropulsive  problems  has  been  demonstrated  using  CRAFT  CFD  flow  solver.  EASM 
extends  baseline  k-c  framework,  validated  to  yield  correct  mean  flow  mixing,  for  a  variety  of  free  shear  flows,  to 
account  for  non-uniform  distributions  in  the  Reynolds  normal  stresses  and  improved  turbulence  realizability  for 
complex  3-D  flows.  Better  comparisons  are  obtained  for  base  pressure  using  a  compressibility  correction  but  at  the 
expense  of  longer  than  experimentally  observed  reattachment  lengths  for  base  flows.  Predicted  streamwise  velocity 
distributions  within  the  recirculation  zone  are  consistently  higher  than  the  measurements  of  Herrin  and  Dutton.  For 
the  lateral  jet  problem,  improvements  using  compressibility  correction  were  significant  and  yielded  a  separation 
zone  size  consistent  with  data.  EASM/J  has  improved  ability  to  resolve  turbulent  flow  structures  in  complex 
vortical  flows,  compared  to  the  k-c  model,  as  demonstrated  by  the  lobes  mixing  study.  The  adequacy  of  the 
compressibility  correction  as  implemented  in  CRAFT  CFD  to  represent  the  flow  features  of  high  speed, 
compressible  flows  still  needs  further  investigation.  This  requires  improved  modeling  of  the  pressure  dilatation  term 
in  the  turbulent  kinetic  energy  equation  of  the  baseline  k-c  model,  and  this  would  further  improve  the  predictive 
capability  of  EASM/J.  Extensions  to  the  c-equation  for  the  problems  described  such  as  curvature  correction  could 
also  improve  comparisons  with  data. 
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